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^— s Abstract 

O^ ' We show that the GPE with cubic nonlinearity, as a model to describe the one dimensional Bose-Einstein 

condensates loaded into a harmonically confined optical lattice, presents a set of ground states which is orbitally 
stable for any value of the self-interaction (attractive and repulsive) parameter and laser intensity. We also 
derive a new formalism which gives explicit expressions for the minimum energy £",„;„ and the associated chemical 
potential /ig. Based on these formulas, we generalize the variational method to obtain approximate solutions, at 
any order of approximation, for £",„!„, /io and the ground state. 

Q^ ' Key words: Bose-Einstein condensates, stability of ground states, analytical approximate formulae, repulsive or 

■^C . attractive interatomic interactions. 
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The Bose-Einstein condensation (BEG) is a fundamental phenomenon connected to superfluidity in liquid 
helium [1]. Its achievement in practice leads to the first unambiguous manifestation of the existence of 
a macroscopic quantum state in a many-particle system. Although the condensates of boson particles 
,^ ■ was firstly predicted by Einstein [2,3] in 1924, BEGs were experimentally realized only in 1995 [4,5]. 

f^ , For these reasons, this phenomenon has attracted the attention of many scientists, in particular, during 

t*^ ' recent years. Nowadays, one of the most interesting problems in cold matter physics is the study of the 

^N I BEG in a potential trap loaded in a periodical optical lattice (for a detailed discussion see Ref. [6] ) . This 

t"*- ■ problem can be described by means of the Gross-Pitaevskii equation (GPE) [7,8,9], for which the BEG 

^D I are characterized by its ground state solutions. 

The ground state solutions of the GPE in an external potential are not necessarily stable, i.e., small 
initial perturbations around a ground state could give rise to solutions collapsing in finite time. In order 
to understand the properties of the BEG, it is important to know how the small-amplitude excitations 
^ . of the ground state evolve in time. One step in this direction was given by Zhang [10], who proved the 

^ ' stability of these solutions in the case of a harmonic potential for a nonlinear attractive interaction. 

Assuming that the harmonic trapping potential has a strong anisotropy (of "cigar-shaped" type), 
the ID limit of the GPE with cubic nonlinearity can be considered as a model to describe the condensate, 
more precisely, by the equation (see [11]) 

^' ^'* + imwVvl/ + Aiz3|vI/|2m/ - Vl cos^ (^x) * = /ip*, (1.1) 



2m dx'^ 2 '"" ' " V f^ , 

where cj > is the oscilator trap frequency, m > is the atomic mass, Vl > is the laser intensity, 
d > is the wavelenght of the laser, /xq G R is the chemical potential and Xuj g M is the self-interaction 
parameter {Xid < when interatomic forces are attractive and Xijj > for repulsive interatomic forces). 
In its dimensionless form the equation (1.1) can be written as 






+ rV^ + AIV-I^ - Vo cos^(aC)V' = mV', (1-2) 
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where, for I := y^n/rnuj, we set ^ = x/l, tp{^) := \/I^(x) and 

2^0 , 2Aiz5 2Vl 2Trl 

raw /aicli Tiuj a 

The solutions of (1.2) can be viewed as standing waves of the time dependent GPE, namely, 

(111 (J u 

V = -7^ + ?'" + ^l"l'" - Vo cos^K)^, (1.3) 

where r := a;t/2 and t is the time. By standing waves we mean time periodic solutions of the form 

M(T,0:=e-'^-^(0. (1.4) 

In this paper we prove the ground state stability (beyond the Bogoliubov approximation) concerning 
the solutions of (1.3), in both cases, attractive and repulsive, and for any value of Vq- This result, together 
with some qualitative properties of the ground states and the general behavior of the minimal energy 
as function of A, completes the essential of Section 2. Furthermore, taking into account the importance 
from the physical standpoint to have explicit formulae for the minimal energy, the corresponding chemical 
potential and the ground state, we develop in Section 3 a new formalism to obtain explicit approximate 
expressions for the above mentioned magnitudes. The present formalism leads to a new and more general 
variational approach. To verify the validity of the method, we compare our solutions with those numerical 
results reported in [12, 13]. 

2. Existence and stability of ground states 

Although the solutions of (1.3) are in general complex valued functions, we can restrict our analysis 
of the existence of ground states only for real valued ones, as we can see by the following lemma, where 
H^(M.) denotes the usual Sobolev space. 

Lemma 2.1: If tp E H^{R) is a complex solution of (1.2), then there exists a real function C/(^) which 
is a solution of (1.2) and a real number 6 such that ip{£,) = e'^C/(^). 

Proof: Assume that ip is a^ complex solution of (1.2) and consider its real and imaginary parts, i.e., 
Tp = u + iv, u ^ 0. Then, it follows that 



d^ 



2 + ^ u + \{u'' + v'^)u - Vq cos (a^) u = ^u. 



(2.1) 
d'^v 



- ^ + C^w + A(m^ + u^)u - Vb cos^ [aC) v^ ^iv. 
dc, 

Now, multiplying the first equation in the above system by v, the second by u and subtracting, we get 



d^v d^u dv du 
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du 
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for some C G R. Since the solutions of (1.2) tends to zero as ^ -> ±00 (see the next Theorem), we 
conclude that C = and hence 

di^ dii d / 7; \ 

0. 

Therefore, v — 7U, 7 7^ and each one of the equations of (2.1) reduces to 

d'^u 
— 777 + ?^« + -^(1 + 7^)m^ - Va cos^ (a^) u = ^lu. 
at, 



Now, considering U{S,) := (1 + 7^)^^^w(C)j it follows that U{^) is a real valued solution of (1.2) and 

where 6 := arctan7. D 

An important property of the solutions of (1.2) is their asymptotic decay at infinity, as asserted in 
the following result. 

Theorem 2.2: Let fx, A, Vb G R be given. If ip e iJ^(R) is a solution of (1.2), then ip e C^{R) and for 
6 e (0, 1) there exists C{6) > such that 

VeeR, MO\<C{S)exp[^{l-6)e/2]. (2.2) 

Moreover, if \ > and /i < 1 — \Vo\, the above inequahty holds for 6 — 0. 

Proof: We proceed as in [14]. Since ip e H^{M.), we have that ip{£,) is a continuous function satisfying 

hm ^(0=0 (2.3) 

ICI-S-+00 

and it follows directly from the equation that -ip" G C(M). In order to prove the exponential decay of ip, 
let 

a(C) ■■=e-^^ + M^m^ - Vo cos2«). 

By Kato's inequality, if z(^) := 1(^5(^)1, we have z" > sign{ip)(p" in the sense of distributions. Therefore, 
—z" + a{£^)z < in the same sense. On the other hand, if we set '(/'o(C) •= C'exp[~(l — (5)^^/2] for 
< (5 < 1 and C > 0, a simple calculation gives 

-V-^,' + a(e)^o - [-5(2 - 5)e -fi + l-S + A|^(e)P - V^ cos^^)] ^o- 

If A > and /z < 1 — \Vo\ set 6 — 0, otherwise assume that < S < 1. Then, for _R > large enough, it 
follows from (2.3) that a(^) > 1 and 

(5(2 - 5)e' -11 + 1-6 + XMOf - Vo cos2(a^) > 0, V|e| > R, 

so that —z" + az < — V'o + '^V'o for \£,\ > R- Moreover, if we choose C > such that z{±R) < ipo{±R), 
then from the maximum principle we infer that 

MO\^m<MO, m>R, 

which implies the exponential decay oi (p, as asserted in (2.2). D 

• Existence of ground states 

We introduce the variational problem which allows to prove the existence and stability of ground 
states for Eq. (1.2). Let 

X := {V G H\M.) ; J (\^\0f + e'lV'(Ol') dx < +c^}. (2.4) 

X is a real Hilbert space if endowed with the following usual inner product 

{<jy\^)x :=y (^'(c)0'(o+cV(c)0(O) ^e 

Then, the associated norm is given by 

3 



Now we define the "energy" E : X ^^ M. and the "charge" Q : X ^^ M. respectively by 

JR JR ^ JR Jm. /r, r\ 

f (2-5) 

and the manifold 



Si:={^eX;Q(^) = l}. 



With these ingredients we look for solutions -0 of Eq. (1.2) that minimizes the energy E among all 
functions in Ei. More precisely, we look for ■0„i„ e Ei such that 

^(V-^.J = min{^(?A); ?Ae El}. (2.6) 

Remark 2.3: Before proceeding to prove that there exist solutions of the variational problem (2.6), we 
shall remember the following well known facts. Let (/jq : R — >■ R be the function 

(^o(e):=^cxp(-eV2). (2.7) 

It is easy to see that ip^ G Ei and that — 'f'o (C) +'?^'/'o(C) — f^iC)- This means that ip^ is an eigenfunction 
of the operator L = —-jfi +C^ corresponding to the eigenvalue Aq = 1. In fact, L has an infinite sequence 
of eigenvalues Aq < Ai < • • •, where A„ = {2n + 1), [n = 0, 1, . . .), and the Hcrmite functions are the 
corresponding eigenfunctions. It is also known that Aq has the following variational characterization. 



Ao = inf{ / (l^'(e)l' + emi)?) dC; V' e El} 



and we can easily verify that Ag = 1 and that the above infimum is actually a minimum attained at ip^^ 
i.e., Ao = llv'ollx == 1- 

We are now in position to prove the existence of ground states of (1.2). 
Theorem 2.4: Let A, Vo G M he given. Then, there exists V'min G Ei such that 

£;(V'„„) = min{s(?A) ; ?A e El}. (2.8) 

Proof: We divide the proof in three steps. 

Step 1: The energy E is bounded from bellow on Ei: 

From Remark 2.3 we know that WipWx ^ 1 for all V" ^ ^i- Hence, if A > 0, we have 

Eiib)>Ufx-\Vo\>l-\Vo\, V^eEi. (2.9) 

On the other hand, from Gagliardo-Nirenberg inequality, there exists a constant Cgn > such that 

||^||^<Cgn||^||^||^'||2, "^yjeH^R), (2.10) 

where || • ||4 and || • ||2 are the standard norms of the spaces L''(M) and L'^{R), respectively. Hence, in the 
attractive case (A < 0), we have for any -0 G Ei, 

\C \C 

EW > M'x + -^Wh -\Vo\> M'x + ^Mx- \Vol (2.11) 

from which we get 

EW> TT^-I^ol, V^eEi. (2.12) 

Id 
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step 2: The variational problem (2.8) has a solution: 

Let E^i„ := mi{E{ilj) ; V € Si}. From Step 1, it follows that i?„i„ G R and from the definition of 
infimum, we conclude that there exists a sequence of minimizing functions {fAnlneN in the manifold Ei, 
i.e., 

Vn e N, BV'n e Si such that hm E{tpn) = i?„i„. 

Assuming that A > 0, we obtain easily from (2.9) that {ipn} is a bounded sequence in X. On the other 
hand, if A < 0, from (2.11) and the Young's inequality, we have 

ll^nll^ < i?(V'n) - ^llV-nlU + l^ol < E{^^) + ^Ujj, + ^^ + \Vo\. 

So, we obtain 

lUnrX<mn) + ^^ + \Vo\, 

from which we conclude that {'tpn} is a bounded sequence in X. Therefore, in both cases, it follows from 
the Banach-Alaoglu Theorem, that there exists a subsequence of ipn that converges to some ■i/'min in the 
weak topology of X, i.e., ipn^^ — ^ V'min- To simplify the notation, we still write ipn for this subsequence. 
Since the embedding X C ^^(R) (see [14,15]) is compact for all 2 < p < +oo, we have 

hm f \MOM^ f H^UOM, 

hm [ \MO\'d^^ [ l^P^UOl'dL (2.13) 

hm / cos^ «) \MO\'d^ = I cos^ «) |^„.„(C)l'de. 



By hypothesis, ■0„ G Si for all n G N, and the second limit above implies that ■0„i„ G Si. Moreover, as 
the norm || \\x is semi-continuous for the weak topology of X , we have 

||^||2,<liminf||V'„||l:. (2-14) 

From (2.13) and (2.14), we conclude that E{ip^i^) — E^,^ and hence V'min is a solution of (2.8). 

Step 3: The function ■0„i„ is a solution of (1.2). 

This follows directly from the fact that E{ip) and Q{ip) are differentiable functionals in X. Indeed, 
from the Lagrange Theorem, there exists /i G M (a Lagrange multiplier) such that 

where E'{ip) and Q'(V') are the Frechet derivatives of E and Q at -0, respectively. Note that this last 
equation is the same as (1.2). This completes the proof. D 

If we denote by Q the set of ground states of (1.2), i.e., 

5:-{V^eSi;i?(7A)=i?_}, 

it follows easily from (2.9) and (2.11) that Q is a. bounded set of X and we have the following properties: 
Theorem 2.5: Let A, Vq G M. 

a) If X > 0, there exists a unique positive symmetric function V'min G Si such that 

^ = {e'V„.„;eGM}. 

b) IfVo = 0, then there exists a positive symmetric function ipmin G G such that ^ i— > V'min(0 ^'*' decreasing 
in the interval ^ > 0. In particular, '!/'mi„(0) = max{V'min(C) ! ^ '= ^}- 



Proof: To prove (a) we proceed as in [16]. Suppose that there are two real functions i/'o, V'l € ^7 V'o 7^ V'l- 
If \tpo\ 7^ \^pi\, define V'l^: = Wi^i + (1 - '^)''PoV^^, where < t^ < 1. It follows that ipi, e Si and 

JR JR /r, 1 r^ 

cos2(aO|V'.(C)l' rf^ = ^ / cos\a()\MOf di + {l-v) ( cos2K)|^i(0I' ^C- 



Since x i~> Ixl** is strictly convex, we have 

\MO\'d^<v I \MO\^d^ + {l-u) I \MO\'dt (2.16) 



Moreover, by differentiating both sides of V'^ — f^V'o + (1 ^ ^jV'i ^-nd using the Cauchy-Schwarz incquahty 
in M^ [ah + cd < \fc^^T^\fWT^), we get 

V-.V^ = i^VoV'; + (1 - ^)V'iV'; < V'.y'KV^^i^^ir^^^KV^)', 

from which it fohows that 

WAi)?<yW^{i)?^{i-v)WAi)?- (2.17) 

As we are assuming that A > 0, it follows from (2.15)- (2.17) that -E(V'i/) < vE^^^ + (1 - v)E^;^ = £"„!„, 
which is impossible by the definition of i?„i„. 

Since |V'i| € Si for i = 0,1, it follows from Kato's inequality that i?(|V'i|) < E{ipi) = £^„i„. So, 
IV'ol) IV"!! € 5 and we conclude by the previous arguments that ji/iol = |^i|- By assuming that V'o(Co) = 
V'i(Co) = for some ^0 G R and taking into account that \tlji\ E C^ (see Theorem 2.2), it follows that 
''/'o(Co) = V'i(Co) = 0, which implies from the uniqueness of solutions of ODEs that ijji = 0. This is in 
contradiction from the fact that ipi G Ei. Therefore, we can assume that ipi — — fAoj with "00(0 > ^or 
all ^ G R. Moreover, since ip{(_): ~ '0(— C) belongs to Q for all ip G G, the unique positive ground state ipo 
is necessarily symmetric. 

To prove (b), let "0 € 5 be a real function and consider V'*(C) the symmetric-decreasing rearrangement 
of |V'('C)I- It is well known (see [9,17]) that ■;/'» is positive, symmetric, decreasing in [0, +00) and 



IMCWd^^ / mOl^d^, l<P<oo, 

-^^ (2.18) 

mo\'d^< / wm'd^. 



Hence, from the first equality of (2.18) with p = 2 we get 1/;* G Si. On the other hand, we have for any 
c> (see [9,17]), 

''{c-eymo\'d^< l{c-e)^\Mi?\di, 



which gives 

c / ' (iV'(e)p - \Mm di< r emo? - \Mm dc (2.19) 

Since the symmetric-decreasing rearrangement is order-preserving, it follows that ip^, also satisfies (2.2). 
Therefore, using the L'Hospital rule we get, for /(^): — \ip{0\'^ ^ l''/'*(OPi 

lim cf /(e)de- lim c^[f{-^c)-f{V~c)]=0 
and consequently 

'eiMoi'd^i^ [emoi'd^. (2.20) 



The first equality of (2.18) with p = 4, together with the second inequality in (2.18) and (2.20) imply that 
E{ip^) = £"„!„, which means that V'* G G- Moreover, for p = +oo, it follows that V'*(0) = max{|-0(^)| ; ^ e 
K}, and the proof is complete. D 

Remark 2.6: Theorem 2.4 asserts that Q is not empty. In fact, Q has infinitely many elements, because 
ii ip e G, then e^^tp E Q, for all e M. Nevertheless, we do not know if Q has only one real positive 
valued function in the case A < 0. However, if there are multiple real valued ground states ip in G, one 
should be aware that the Lagrange multiplier /x might depend also on -0. Indeed, multiplying both sides 
of Eq. (1.2) by ip, we get 

^^W^E^^„ + ^ [lijiOM. (2.21) 



As we are going to see in the study of stability, it is important to notice that the set of Lagrange 
multipliers is bounded. This is immediate because 

Hn < \E^..\ + ^IIV^Ill < \E^..\ + ^-^Mx 

and because G is bounded in X. 

Remark 2.7: Theorem 2.4 states that Eq. (1.2) admits ground state solutions in both cases: attractive 
(A < 0) and repulsive (A > 0). In the repulsive case, Eq. (1.2) presents a different behavior when 
compared with the classical NLS. In fact, assume that (p € H^{R) is a nontrivial solution of 

-ifi" + X\ipf(p^ Hip. (2.22) 

If we multiply both sides of (2.22) by </?(^) (respectively — ^</?'(C)) smd integrate on K, we get 



Hence, 



\^'iO\'dC+^ flpiOl'd^^O, 
^ Jr 

which is impossible if A > and f ^ 0. Therefore, (2.22) has no nontrivial solution in H^(M.) if A > 0. 
• Stability of ground states 

In order to prove the stability of ground states of (1.2), let us consider the Cauchy problem 

i^ = ^E'iv), ^;(0,C)=«o(0, (2.23) 

where E' is the Frechet derivative of E in X, i.e., 

^E'{v) :- -0 + ev + X\v\^v - Vo cos^ «) v. 

It is well known [18, 19] that (2.23) has a unique solution that is global in time, i.e., for any vq G X, 
there exists a unique v € C([0, +00), X) satisfying (2.23). In particular, ii 4' ^ Gi the unique solution u 
of (2.23) such that u(0,^) — V'(0 is the standing wave given by 

^i(T,0 = e^'^"^(0- 

Saying that ip is stable means that if the initial datum v^ of Eq. (2.23) is close enough to ip, then the 
trajectories w(t, •) remain close to the set Q, a.s t varies in R. More precisely. 



Definition 2.8: Wc will say that Q is stable if, for each £ > 0, there exists S > such that if vq E X 
satisfies ini^^g \\vq — ip\\x < S, then the solution of (2.23) satisfies 

sup inf \\v{t, •) - e^^'^^tpWx < e. 

Before proving the stability oi Q, we state and prove the fohowing lemma about the compactness of 
the set Q, which will be needed in the proof of the stability result. 

Lemma 2.9: The set Q is compact and weakly sequentially closed in X. More precisely, if {tpnin&n is 
a sequence ofQ, then there exists ip & G 3,nd a subsequence {ipn^}keN such that "tpn^ — > V' i^ ^- Also, if 
{'>pn}neN is a sequence ofQ and ipn ^" i^ ™ X-weakly, then ■)/'„-> -0 strongly in X. 

Proof: Let {tpn}n<£N be a sequence of Q. Since Q is bounded in X , there exists a subsequence {'0„j.}fcgN 
and ip E X such that ipn^. -^ ip m X-weak. By the compactness of the embedding X C LP{M.) for 
2 < p < oo, we conclude that ip £ J^i. Now, arguing as in the step 2 of the proof of Theorem 2.4, we 
obtain that ^ E Q. To conclude the proof, it remains to show that ?/;„^, — > -0 in X strongly. To see this, 
it is enough to remark that ||'0nfc||x -^ llV'llx, but this follows from (2.13) and the fact that 

UnAW = ^».n -^ I |^n.(Ol'de + VJ COS^ «) |^„,(C)PdC. (2.24) 

^ JR JR 

This proves that Q is compact and the remaining claims follow easily. D 

In order to prove the stability of Q, we need the following well known conservation laws that hold 
for all solutions of (2.23): Assume that v{t, ^) is a solution of (2.23) such that v{0, ^) ~ fo(C)- Then, for 
each r e M, we have 

Q{v{t, •)) = Q{vo) and S(i;(r, •)) - E{vo). 

Concerning the above identities, the first one is known as the "conservation of charge" and the second 
one is the "conservation of energy". The first one can be obtained multiplying Eq. (2.23) by — iti(r, ^) 
(i.e., the complex conjugate of iv{T, ^)). By the same way, one obtains the conservation of energy, but in 
that case we should multiply the equation by ^. 

We are now in position to prove the stability of Q. 

Theorem 2.10: The set Q is stable in the sense of the Dehnition 2.8 . 

Proof: Arguing by contradiction, we proceed as in Cazenave-Lions [20]. If G is not stable, there exists 
eo > such that for all integer n e N, we can find vq^ G X satisfying 

r„ := inf ||wo« - ip\\x < - (2.25) 

and 

sup inf ||«„(T,.)-e-'^W-^||x>£o, (2.26) 

where f„ G C(M; X) is the unique solution of (2.23) with initial datum won- 

Let 4'n & G such that r„ < ||i;on — ipnWx < Tn + l/n. Since G is bounded in X and the sequence 
IJ.{ipn) is bounded in M (see Remark 2.6), there exists {ipocfJ-ca) & X xM. and a subsequence still denoted 
by (0„,/z(i/'n)) such that 0„ converges to tpoo weakly in X and strongly in L*(M.) n L^(M), while /i(V'n) 
converges to /ioo in M. By Lemma 2.9 we know that 0„ — > i/'oo in X and tpoo G G- Also, from the fact 
that ll^on — V'nilx — > 0, we may infer that won converges to -000 in X and in L*{M.) n _L^(M) as well. 

On the other hand, one may observe from (2.26) that there exists t„ G M such that 

inf |1«„(t„, •) - e-'^'W^-^IU > ho. (2.27) 



Setting V'n := e'^('^")'^"w„(T„, •), it follows from the conservation of charge and energy that 

Qi^n) = Q{vOn) -^ Qii^oo) = 1, ^(V'n) = E{vo„) ^ E{^oo) = ^„i„ aS 71^ CO. 

All this means that {V'nIraeN is a bounded sequence in X and there exists -000 G -^ and a subsequence 
(still denoted by {V'nInGN) converging to ipoo weakly in X and strongly in L'*(R) n L^(IR). Therefore, 
V'oo € Q and E{tpn) — > ^(V'oo)- Again, invoking relation (2.24) we observe that i/'n ^ V'oo strongly in X 
and this is a contradiction with (2.27). This finishes the proof. D 

• The minimal energy as function of A 

In order to explicit the dependence of the minimal energy relatively to the parameter A, let us denote 

i;„.„(A):=min{£;A(0);0eSi}, (2.28) 

where Ex is the energy functional introduced in (2.5), and Qx := <ip (^ Y,i ; Ex{ip) — ii'„i„(A) > the set of 

corresponding ground-states. Theorem 2.4 assures that i?mi„ is well-defined as function of A G R and it 
is easy to see that it is strictly increasing. Indeed, for /i > and tp G Gx+h we have 

£^„.„(A) < Exii^) - ^„„(A + h)- f^mi < ^_(A + h). (2.29) 

Proposition 2.11: _E,„i„(A) is a strictly increasing and concave function such that 

lim .B„,„(A) =±cx). (2.30) 

A— >±oo 

The proof of the Proposition relies on the following: 

Lemma 2.12: For each a,b gM., a < b, the set Ua<A<6 ^^ ^'® bounded in X. More precisely there exists 
a constant Ca.b > such that 

Mx<Ca,b V^e U Gx. 

a<X<b 

Proof: It suffices to prove for a < < 6. Let A G [a, b] and ^p £ Gx- Then 

E^Ub) > E„U\) = ExW > MJc + ^\mt - \Vo\- 
If A > 0, then HV'llx < E^i„{b) + |Vb|. If A < 0, it follows from Gagliardo-Nircnberg inequality (2.10) that 

2 , •^^'gn,, II I I ^ II ,||2 . "C'gn, 



and we have 



E^Ub) > Mx + -^Mx - \Vo\ > Urx + -^Mx - \Vo\ 



Ufx < 2 K.„(6) + ^$^ + \Vo\] ■ a 



Proof of Proposition 2.11: It follows from (2.29) that 

£^„.„(A) < liminf £;„„„(A + h). (2.31) 

Let ^ eGx- Then, 

^„.„(A + h)< Ex+h{'4>) = E^JX) + ^\mi yh e R, (2.32) 



from which we obtain 

limsupi;„,„(A + h) < S„,„(A). (2.33) 

From (2.31) and (2.33) we eonclude that 

hm E^,JX + h) = -B„,„(A). (2.34) 

On the other hand, ii h < and ip G Gx+h, we have 

^„.„(A) < Ex{ib) - £^„.„(A + h)- ^ll^llt (2.35) 

from which, using the Gaghardo-Nirenberg inequahty we obtain, 

hC 

E^U^) < i?_(A + h)- ^II^IU. (2.36) 

By choosing and fixing a G M such that a < X + h < X, it follows from Lemma 2.12 that there exists a 
constant C'a.x such that 

hC 

i?„.„(A) < E„UX + h)- -^Ca,x, 



from which we get 

So, (2.33) and (2.37) give 



liminf S„,„(A + h)> S„i„(A). (2.37) 



lim E^,JX + h) =S,„„(A) 



and we conclude from (2.34) that i?„i„ is a continuous function. Moreover, from (2.32) we have 

-E„,„(A + h) + .B„,„(A -h)- 2£;.„„(A) < VA, /i G K, 

which, under continuity, is sufficient to assure the concavity of i?„i„. 

In order to prove (2.30), consider ipk G Gk, fc G N. Then, using (2.29) and (2.32) with ft, = 1 we have 

lUk+iWt < E^Jk+1) - E„,Uk) < lUkWi ykeN 
from which wc obtain for all n G N 

-. n+l -. n 

fe=l fe=0 

Now, arguing by contradiction, assume that there exists a positive constant C such that i?„i„(A) < C 
for all A G M. Then, it follows from the first inequality in the above expression that the sequence 
{V'nin converges to zero in L*{M.). But, repeating the arguments used in the first step of the proof of 
Theorem 2.4, we conclude that the sequence {tpk}k is bounded in X, so that, passing to a subsequence 
if necessary, we have that ipk converges to ip weakly in X. Since the embedding X C L^(R) is compact 
for 2 < p < oo, we conclude, by choosing p = 4 that t/j ~ 0, and by choosing p ~ 2 that ip G Si, which is 
absurd. Since the same arguments apply for ip-k G G-k, fc G N, we finish the proof. D 

Remark 2.13: Concerning the dependence of the chemical potential relatively to A, it follows from the 
formula (2.21) and the characterization oiGx given by Theorem 2.5 that we can also define, for A > 0, the 
function /i„i„(A). However, as reported in Remark 2.6, it is unclear that /^„i„(A) be uniquely determined 
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if A < 0. No matter whether or not it be uniquely determined, the fact that X(^ii\{ip) — _E„i„(A)) > for 
all A G M and for all ip E Gx, allows to state that 

lim ^i\{ip\) = ±00, yxpx e Gx. 

A— ^■ioo 

Corollary 2.14: Let ipi e G\,, i = 1,2. If Ai < A2, then \\tlj2\\i < ll^/'ilU- 
Proof: Let h ^ X2 - Xi- Then, from (2.29) and (2.32), it follows that 

1|, , |,4^ i^n,.n(A2)-g„„,(Al) ^ 1„, „4 

2llV'2|i4< j^ < ^imu- □ 

3. A new method to obtain approximations for the minimal energy 

The results presented in the previous section, although mathematically rigorous, do not provide 
sufficiently precise quantitative information for some relevant physical quantities. Since it seems to be 
not possible to calculate exact explicit solutions of Eq. (1.2) (in particular, the ground state solution V'min), 
the exact values of such quantities cannot be expressed in terms of the known parameters. Therefore, 
it would be useful to obtain some kind of explicit formulae through which one could approximate them. 
Indeed, these are outcomes mainly interesting from the physical point of view. In this sense, there exist 
some well known methods that were already applied to the problem we are dealing with (see [12,13,21]). 
Nevertheless, as we shall see below, a new interesting and powerful approach can be developed. 

Let {(^a}agr a family of functions in Ei such that Xi-^ (px defines a differentiable curve in X. Then, 

MO-^^xiOd^^O VAeM. (3.1) 

Since the energy functional E\ is differentiable in X, it follows from the chain rule that A h- > E\{ip\) is 
also a differentiable function and 

^Ex{vx) - {e'x{vx) I ^^a) + \j^ Wxm'di. (3.2) 

where (-j •) denotes the duality product between X and its dual X* . 

From now on we assume that, for any A G R, one can choose a real ground state ip\ S G\ such that 
A I— >■ i/^A defines a differentiable curve in X. In this case, 

i?„,.n(A):= I Wx{0?d^+ I e\MO?d^+^ f\MO\'d^-Vo f cos\aO\MO\'d^ 

JR JR ^ JR JR 

is differentiable as a function of A and, as E'^{ip\) — iJ,ifi\, we get from (3.1) and (3.2) 



|j£.,„(A)^l/j,yoi«d?.i|lv, 



AII4 



(3.3) 



and we have the formula 



i?_(A)=i?_(0) + iy UsWtds. (3.4) 



On the other hand, if we denote /i„i„(A): — iJ,{ip\), it follows from (2.21) that 



^,,.ux)-m\t + ^,j-,m\i (3.5) 



and we get by integration on A the formula 



M„.n(A) = Ai_(0) + MA||^,||4 + y UsWtdA. (3.6) 



Notice that £'„i„(0) — /i„i„(0) and, in the case Vq = 0, we have from Remark 2.3 that £^„i„(0) 
M™(0) = 1. 
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Remark 3.1: Formulae (3.4) and (3.6) express the minimal energy E^i„ and the corresponding chemical 
potential y^n,i„ as functions of A, which also depend explicitly on the L'^-norm of unknown ground states. 
However, by eliminating this explicit dependence, we can obtain an exact formula relating these two 
quantities. More precisely, by using (3.3), we can rewrite (3.5) as 

^...(A)^2A^^,„(,) + ,|^^^^„(,), 



from which we get easily 



i i5„„(A)-i/„.„.W<i.Uo. 



Hence, there exists a constant C such that 



i?„.„(A) 



M.„.„(s) ds == C, VAe 



Since i<^„i„(0) = /imin(O), it follows that C = and we have the identity 



^„.„(A) = - / Mn.„(s)ds, VAe 
A Jo 



(3.7) 



The formula (3.4) and (3.6) can be used to obtain explicit approximate functions depending on A 
for the minimal energy i?„i„ and the corresponding chemical potential /i„i„, by choosing appropriate trial 
functions that generate differentiable curves in Si. 

Motivated by inequality (2.2), we consider the trial functions (p^ G ^i defined as 



MO--^ 



■ exp(— K^ 



A direct calculation shows that the function k h- > E\{(J)k) is given by 



^^("^'' + s + ^-fi'—' 



a 
2k, 



We can show that, for all A G R, the value of k that minimizes the function k m- E\{(J>k) is given by the 
largest solution (in fact, unique solution if o^Vq < 1) k(A) of the following transcendental equation 



^3/2 f «:i/2 



A 



4V^ 



Voa^ 



■ exp 



2k 



(3.8) 



In order to obtain explicit approximate formula for the minimal energy E^^^ (A) and the correspond- 
ing chemical potential /i„i„(A), we introduce the functions 



A*app(A): - M™(0) + I UyxWt + ^ 11^.11^ rf.s 



where, in this case, 



v\{0-^ 0k(a)(O 



2k(A) 



1/4 



exp(-«:(A)C'). 
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By a straightforward calculation we get 



SV"" Jo , , 

\ (3.9) 



M.pp(A)=A*™(0) + ^|AvMA)+y^ ^W)ds 



Moreover, by arguing as in Remark 3.1, we can show that 



1 '•^ 



£^app(A) = ^ / fJ-.pAs)ds. 

From the above identity we can easily relate the Taylor coefficients _E„ of -E^pp with the ones of /i^pp- In 
fact, we have Hn = (n + l)En for all n G N. 

A relatively simple situation is the one at which the optical lattice potential is not present (Vq = 0). 
In this case, (3.8) has a unique solution and if we define cr(A): = ■\/k(A), the equation (3.8) with Vb = 
can be written as 

40r 4 ^ ^ 

By differentiating this equation implicitly with respect to A, we get easily the following properties of the 
function a{X): 

Lemma 3.2: The function a : M. ^ R is C°° , positive, strictly decreasing, convex and satisfies the 
foUowing properties: 

cr(0) = \/2/2; lim cr(A) == +oo and lim cr(A) = 0. 

A— ^ — CXD A— ^ + oo 



A , .... fV^-'^' 



More precisely, 

ct(A) ^ '-^ as A ^- — oo and cr(A) ^ ( -^^ ) as A ^ +oo. D 

Ay'TT V A / 

As consequence of the previous lemma, we can show that _E^pp(A) satisfies the general properties of 
E^in{X) as reported in Proposition 2.11. More precisely. 

Corollary 3.3: The functions E.^pp{X) and /i^pp(A) defined in (3.9) are C°° , strictly increasing and 
concave. Moreover 

{ i?app(A) and a*bpp(A) arc O(-A^) as A -^ -oo, 

{ ^ 

\ Eapp(A) and A*app(A) are 0{\^''^) as A ^ +oo. 

By implicit differentiation the Eq. (3.10) on A, we get 



^^(0) = ^, a'(0) = -— ^, a"(0) = -— ^, a"'(0) = --— ^, aW(0) 



2 ' ' ' 16V^' ' ' 128^^2' 5127rV^' ' ' 16384^2^' 

(3.11) 
Then, for A small enough, we can consider the approximation 

,,, V2 A 3A2 A-^ 45A4 

o-(A) ~ ^r ' 



2 16V^ 2567ry2 10247rV^ 3932167r2\/2' 
By substituting the above approximation in (3.9) we have, for A small enough, 

A A2 A3 A4 9A5 



£'app(A) Si 1 



A*app(A) W 1 + 



2V2^ 647r 5127^^/2^ 81927r2 7864327r2\/2^' 
A 3A2 4A3 5A4 54A5 



(3.12) 



27r 647r 5127rV2^ 81927r2 7864327r2^/2^■ 
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Remark 3.4: The above results are valid for attractive (A < 0) as well as for repulsive (A > 0) interatomic 
interaction strengths. By using (3.3), (3.5) and (2.7) we can show that E'^.^{0) = E'^ (0) = l/2-\/27r, 
which implies that the formulffi given by (3.12) coincide with ii^„i„(A) and /imin(A), respectively, up to 
first order terms in A. Notice that, up to second order terms in A, the approximate chemical potential 
/iapp can be written as 

/x.,,(A)«l + ^-eA2, (3.13) 



where e = jj- ss 0.0149207. 

64-7r 

It is noteworthy to compare Eq. (3.13) with the one obtained in [12] by a perturbative method. 
Considering that our dimensionless parameters are in fact twice the ones used there, both formulae (3.13) 
and Eq. (31) in [12] coincide, except for the values oi e: e = 0.0149207 and e = 0.016553 respectively, 
which are also very close to each other. Moreover, Fig. 1 displays a comparison between the values of 
Mapp from Eq. (3.12) (solid line) and those obtained by perturbation theory (dashed line) [12]. Also and 
for sake of comparison, the numerical evaluation of Eq. (1.2) is shown by full stars. From the figure is 
observed that Eq. (3.12) increases the accuracy of the solution with respect to the perturbative method. 
In the scale of the figure no significant differences are observed between the numerical solutions and the 
calculated values using Eq. (3.12) for the interval range jAj < 8. Nevertheless the perturbation method 
gives a large error for jAj > 6. 
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Fig.l: Dimensionless chemical potential /i = 2iiqITvjJ as a function of A = 2\i£i / {Iftio) . Solid 
line: calculation following Eq. (3.12). Dashed line: perturbation theory from Ref. [12]. Stars: 
numerical evaluation of Eq. (1.2). 

Furthermore, in [13] a closed expression for the order parameter is given by (in the case Vq — Q) 



1 ^2 

V'(C) = i;=cxp(- — ) 



A 
2tt ji 



72/2 



exp(-e2(l-z2), 



dz 



In our approach, we propose the function of Ei: 



^A(e)=('^^V^'cxp(-«(A)e^), 



(3.14) 



where n{\) is the unique root of the equation (3.8) with Vq = 0. Recalling that fi;(A) = cr(A)^, it follows 
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from (3.11) that, for A small enough, 



;,(A)«i-lW-A + ^A2, 
^ ' 2 16 V TT 1287r ' 



V^) 



1 



^ 



A + 



^ 



^ 320F 10247rV2 
and (3.14) can be approximated up to second order terms in A by 



A^ 



'/'app(0 



1 f e 

TT^exp 



1 



A2 



V2X 
32V7r ^ 10247r 



exp 



A 



A2 



8V2^ 1287r 



(3.15) 



It is also interesting to notice that, for all A G 

1 A 



A2 1 



> -. 



2 8v^ 1287r - 4' 

v^A A2 1 

1- -^-^H > -, 

32^^ 10247r " 2' 

which implies that Ip^^p is a positive function of X for any A G M. 

Remark 3.5: It follows from the properties stated in Lemma 3.2 and the Labesgue Theorem that 

lim \\(f>x - ^o\\x = lim \\ip\ - ipo\\x = 0, 

A-J-O A-5-0 

where (always assuming that Vq = 0) (po, fx and ifx are given by (2.7), (3.14) and (3.15), respectively. 
Therefore, from Theorem 2.10, for A small enough and up to a change of phase, the unique solutions of 
Eq. (1.3) with initial data ipx and ipx respectively, remain close (in the sense Definition 2.8) to m(t, ^) := 
e^'^</?o(0 in X, for ah time t G R. 

4. Conclusions 

Our main results in the first part of this work concern qualitative properties of the minimal energy 
solutions of ID GP equation with cubic nonlinearity in a harmonically confined periodical potential. We 
prove the existence of ground states for any A and Vq (Theorem 2.4). Regardless the value of the laser 
intensity Vq, such ground states have a Gaussian-like exponential asymptotic behavior, as was pointed 
out in Theorem 2.2. For A > (repulsive interatomic forces), we prove that there exists a unique positive 
and symmetric ground state i/'min, which is decreasing for ^ > if Vq ~ 0, and that any other solution 
of the problem differs from it in a phase factor (Theorem 2.5). We also prove that, independently of 
the value of Vb and for any A G M, the set G\ of ground states is orbitally stable (Theorem 2.10). An 
important consequence of this fact is the stability of those physical magnitudes, which are described by 
operators defined in the Hilbert space X (see (2.4)). This result is particularly related to the superfiuidity 
properties, among other physical phenomena, of the harmonically confined condensates loaded in optical 
lattices [22]. 

In the second part (Section 3) we present a new and simple method to construct formulae (see (3.4) 
and (3.6)) that allows to approximate the minimal energy, the corresponding chemical potential as well 
as the ground state. The functions described by these formula (see (3.9)) preserve some global properties 
of -E„i„(A) and /i„i„(A), as pointed out by Proposition 2.11 and Corollary 3.3. In the case Vo — 0, we 
obtain approximations for ii^„i„(A) and /i„,i„(A) as Taylor polynomials of E^^p(X) and /i^pp(A) respectively 
(see (3.12)). 
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